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ABSTRACT 


This report consists of two parts: 1. A brief description of the implementation 
of a new equation solving and eigenvalue extraction package into BOSOR6 (a program 
for the analysis of hybrid structures of revolution) and 2. a new User's Manual 
for BOSOR4, which has appeared as a chapter in STRUCTURAL MECHANICS SOFTWARE 
SERIES - VOL. 1, edited by Nicholas Perrone, Walter Pilkey and Barbara Pilkey and 
published by the University Press of Virginia in 1977. 


Nov. 30, 1977 


Mr lurcfy 
Dear Mrs.—Taylor, 


Enclosed are the originals for the final report for 
which DDC expressed a need for a better manuscript from which 
to make copies. 


Sincerely yours, 
Gri DAA 


David Bushnell 


Contract NO0014-76-C-0692 
Part I 


Incorporation into BOSOR6 of a New Simultaneous Equation Solver and Eigenvalue 


Extraction Subroutines, 
ABSTRACT 


Subroutines written by Frank Brogan for the STAGS computer program [1] were 
incorporated into BOSOR6 [2] and check cases were run to debug the implementation. 
Most of the effort in this task was devoted to restructuring the BOSOR6 mass 
storage data files to make them compatible with the out-of-core equation solver 


and eigenvalue extraction package. The relevant subroutines are listed with their 


purposes defined on the following two pages. 
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DATA MANAGER 
FMM Purpose: 
SPACE 
GASP 


gw 
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MATRIX DECOMPOSITION ROUTINES 


Purpose: 


GLOBAL MATRIX ASSEMBLY 


Purpose: 


Keeps track of data files stored 
on random access mass storage 
device. 


Special purpose calls to FMM. 


Decomposes a symmetric matrix into 
lower and upper triangular factors. 
The matrix is stored in blocks on 
random access mass storage device. 
The data manager FMM is used to 
store and retrieve blocks of the 
symmetric matrix to be decomposed, 


Assembles the total stiffness matrix 
from the element stiffness file. This 
assembly is accomplished by means of 
a previously calculated skyline vector. 


EIGENVALUE EXTRACTION ROUTINES 


Set up arrays of pointers to addresses 
of data in core or on mass storage 
device. These data include the element 
stiffness blocks, mass blocks, load- 
geometric blocks, assembled stiffness, 
mass, and load-geometric matrices, 
eigenvectors, skyline vector, and sign 
vector. 


Extract eigenvalues for modal vibra- 
tions or bifurcation buckling. 


The CLUSTER program has been designed 
to solve the generalized eigenvalue 
problem KX = ,BX when the equation 
system is very large and the number of 
eigenvalues desired is small, usually 
less than ten. A relatively simple 
combination of simultaneous inverse 
iteration, Chebyshev acceleration, 

and reduced eigenvalue solution pro- 
vides an efficient method for most 
problems of this type. In particular, 
CLUSTER avoids additional shifts of 
the eigenvalue spectrum and consequent 
recomputation of the factored stiff- 
hess matrix except in a few special 
circumstances. 
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REFERENCES 


PART I 


1. STAGS Theory Manual; Currently being written by B. O. Almroth and F. A. Brogan. 


2. Bushnell, D., "Stress, Buckling and Vibration of Hybrid Bodies of Revolution, 
Vol. II: User's Manual for BOSOR6", LMSC-D501504, March 1976. 
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Contract N00014-76-C-0692 


PART II 


New BOSOR4 User's Manual 


ABSTRACT 


The following User's Manual has appeared as Chapter 1, pp. 11-142 in_Structural 


Mechanics Software Series, Vol. 1, edited by N. Perrone, W. Pilkey and B. Pilkey. 
This manual, an update of "Stress, Stability and Vibration of Complex Branched 
Shells of Revolution: Analysis and User's Manual for BOSOR4", LMSC report 
D243605, March 1972, was written, laid out, edited, and typed under this contract 
in response to a request by one of the contract monitors. This task was per- 
formed in lieu of the second objective stated in the July 1976 annual report for 
CONTRACT NQ0014-76-C-0692. The purpose of this second task was to make the 
BOSOR6 program capabilities conformable with the capabilities implied by the 


user documentation. 
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Loads and Temperatures Distributed over the Surface of Segment “ISEG" 


500 Continue 
@ NLTYPE, NPSTAT, NISTAT, NIGRAD 


a 


Definitions of Input Variables 


NLTYPE = control integer for type of loading: 


0 = no pressure, surface traction of temperature distribu- 
tion on this shell segment. 

1 = pressure and/or surface traction, but no temperature 
distribution on this shell segment. 

2 = temperature distribution, but no pressure or surface 


traction on this shell segment. 
3 = pressure and/or surface traction and temperature distri- 


bution on this shell segment, 


number of meridional stations in this segment for which pres- 
sure and surface traction components will be read in. If 
INDIC = 3 or 4 and NLTYPE = 1 or 3 NPSTAT mist be greater 
than or equal to 2,and less than or equal to 20. The NPSTAT 
= 0 option can be used only for INDIC = -2, -1, 0, 1, and 2. 


il} 


NPSTAT 


Sacma (2) NTSTAT = number of meridional stations in this segment for which 
ave f ~~ temperature rise coefficients Tl, T2, and T3 will be read in. 
2 Same discussion applies here as for NPSTAT. 


NIGRAD = control integer for type of thermal gradient through the 
shell wall thickness: 
2 Ss 
l means T(s,z) = Tl(s) + T2(s)*z + T3(s)*z A, 
2 means T(s,z) = Tl(s) + 12(s)*27 36S) Ly 
z 


4 3 means T(s,z) Tl(s) + T2(s)*exp(z*T3(s)) 
where z is measured from the reference surface positive to 
the right of increasing meridional arc length, s. In 
Tables A3 and A4 the function T(s,z) is called "FUNCT." The 
j actual temperature magnitude is given by T(s,z)*TEMP or 
T(s,z)*DTEMP for INDIC = -2, -1, 0, 1, or 2, and for each 
BOSOR4 Data Deck circumferential harmonic by T(s,z)*TDIST(L,ISEG) for INDIC 
equal to 3 or 4. 
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